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ABSTRACT 

A reduced dynamical model is derived which describes the interaction of weak inertia- 
gravity waves with nonlinear vortical motion in the context of rotating shallow-water flow. 
The formal scaling assumptions are (i) that there is a separation in timescales between 
the vortical motion and the inertia-gravity waves, and (ii) that the divergence is weak 
compared to the vorticity. The model is Hamiltonian, and possesses conservation laws 
analogous to those in the shallow-water equations. Unlike the shallow-water equations, 
the energy invariant is quadratic. Nonlinear stability theorems are derived for this system, 
and its linear eigenvalue properties are investigated in the context of some simple basic 
flows. 
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1. Introduction 

An important area of research in geophysical fluid dynamics is the development of 
simplified models. Simplified models are reduced models of the full hydro dynamical equa- 
tions of motion that retain the essential physics of the problem under investigation, yet 
are reduced in complexity and are therefore more amenable to theoretical analysis. 

A generic property of atmospheric and oceanic dynamics is the existence of nonlinear 
vortical motion, evolving on an advective timescale, together with relatively fast acoustic 
and gravity wave motion. While both kinds of motion can be identified observationally, it is 
generally the case that the vortical motion overwhelmingly dominates the fast wave motion 
for both meso-scale and large-scale dynamics. This fact has led to a class of simplified 
models known as balanced models, which filter the fast wave motion and describe only 
the nonlinear vortical motion. Examples of such models include the barotropic vorticity 
equation, the quasi-geostrophic equations, and the semi-geostrophic equations (see e.g. 
Pedlosky 1987). 

Yet while the fast wave motion is observationally found to be relatively weak, it is 
nevertheless ubiquitous. Moreover, there are strong theoretical reasons for believing that 
balanced models cannot be exact (even in principle), and that some degree of fast motion 
is inevitable. There is, therefore, considerable interest in the problem of the interaction be- 
tween (slow) nonlinear vortical motion and (fast) weak waves. Clearly it would be desirable 
to avoid dealing with the full primitive equations of motion in theoretical investigations of 
this nature. Yet simplified models describing this kind of interaction do not exist, except 
for very special cases. 

In this paper we address this problem of constructing simplified models describing the 
interaction of nonlinear vortical motion with weak waves under a separation in timescales. 
A traditional approach is through perturbation methods, employing an asymptotic expan- 
sion with the timescale ratio between the wave and vortical motion (i.e. the Mach, Froude, 
or Rossby number) serving as the small expansion parameter. However, such methods are 
generally carried only to leading order. This is, perhaps, not surprising in light of the 
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fact that the higher-order corrections in traditional expansions will generically exhibit un- 
bounded secular growth in time (Warn et al. 1995). Yet to obtain non-trivial coupling 
between the wave and vortical motion, one has to go beyond leading order. 

A rather different approach is to argue that one should maintain the fundamental 
dynamical properties of the original system, notably its conservation properties. The im- 
portance of maintaining conservation properties in dynamical simplifications has been long 
recognized (Lorenz 1960, Sadourny 1975). More recently, it has been realized (e.g. Salmon 
1983) that the most efficient way to do this is through maintaining the Hamiltonian struc- 
ture of the dynamics: the simplifications in question invariably affect the conservative part 
of the dynamics, which for nearly every model in geophysical fluid dynamics is Hamiltonian 
(see e.g. Benjamin 1984, Shepherd 1990). 

Ideally one would like to be able to satisfy both requirements, Hamiltonian structure 
and formal accuracy. However, so far they have proved to be incompatible beyond leading 
order. Arguments can be mustered for one approach or the other. We do not wish to enter 
that debate here, but simply choose to follow the Hamiltonian approach, guided by scaling 
arguments, in what follows — recognizing its possible limitations. 

The system we shall investigate is rotating shallow- water flow. This is chosen because 
it is arguably the simplest system in geophysical fluid dynamics that contains the essential 
physics of nonlinear vortical motion coupled to fast (in this case inertia-gravity) waves. 
Yet the approach may be expected to have implications for other systems having the same 
essential characteristics, such as vortex/acoustic- wave interactions in compressible flow. 
Farge & Sadourny (1989) proposed a weak- wave model for shallow- water flow, but their 
derivation was rather ad hoc and their system has a number of unsatisfactory features. We 
are aware of no other study along these lines. 

The plan of the paper is as follows. In Section 2 the properties of the shallow-water 
equations are briefly reviewed. Our weak-wave model is derived in Section 3, in a way 
that guarantees its Hamiltonian character. We then describe some of its properties: sym- 
metries and conservation laws (Section 4); pseudoenergy and pseudo momentum invariants 
and nonlinear stability (Section 5); rigorous upper bounds on gravity- wave emission from 
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unstable parallel vortical shear flows (Section 6); linear eigenvalues and instability (Section 
7). Some of the limitations of our approach are discussed in Section 8. 



2. The shallow- water system 

The shallow-water equations, describing the dynamics of a shallow, homogeneous, 
rotating fluid, are given by (e.g. Pedlosky 1987) 

v t + (v ■ V)w + fz x v = -gVh, h t + V • (hv) = 0, (2.1) 

where v(x, y, t) is the (horizontal) velocity, h(x, y, t) the fluid depth, g the (constant) grav- 
itational acceleration, / the (constant) Coriolis parameter, z the unit vertical vector, and 
V = (d/dx,d/dy) = (d x ,d y ). We assume either an unbounded domain with appropri- 
ate decay conditions at infinity, or a doubly-periodic domain, or a combination of the 
two (i.e. unbounded in one direction and periodic in the other). This system has the 
following integral invariants: energy Tt = J ^(h\v\ 2 + gh 2 ) dx; zonal absolute momentum 
M = J h(u—fy) dx, where v = (tt, v); angular momentum J = J h[xv—yu+^(x 2 +y 2 )] dx; 
and a family of Casimirs C = f hC(q)dx where C(-) is an arbitrary function and q = 
[f + z • (V x v)]/h is the potential vorticity. The Casimir invariants are a consequence of 
the material conservation of potential vorticity, q t + v ■ Vq = 0, which follows from (2.1). 
The system (2.1) can be written in the generic Hamiltonian form (e.g. Shepherd 1990) 

a t = J^, (2.2) 

where a = (v, h) T is the state vector of the dynamical variables, H is the Hamiltonian, 
and the cosymplectic form J is given by 



/ q -d x \ 



J = 



\-d X -dy / 



(2.3) 
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If one linearizes the dynamics (2.1) about the rest state v = 0, h = H = const., the 
dispersion relation is third order: one root has frequency ui = 0, corresponding to the vor- 
tical motion; and there is a pair of roots with frequency u = ±/a/1 + gHn 2 / f 2 L 2 , where 
L is a horizontal length scale and k the 0(1) dimensionless wavenumber, corresponding 
to inertia-gravity waves. Allowing weak nonlinearity, the vortical motion has a timescale 
L/V, where V is a typical velocity. The ratio of the inertia-gravity wave timescale to the 
vortical timescale is thus e = V/\oj\L = RB + B 2 , where R = V/fL is the Rossby 
number and B = fL/^/gH. Evidently there is a separation in timescales when e C 1. If 
one further assumes that the timescale of the dynamics is advective (so that the gravity 
waves are filtered), one obtains balanced dynamics: quasi-geostrophic dynamics for R <C 1 
with B arbitrary, and barotropic dynamics for B <C 1 with R = 0(1) (or equivalently 
F<1, where F = V/y/gH is the Froude number). 

It may be noted that the energy of the full shallow- water equations (2.1) is cubic in 
the dependent variables, which is a consequence of the fact that the flow is divergent. (The 
same situation arises with compressible flow.) This presents some serious difficulties with 
respect to theoretical investigations. First, the partitioning of energy between different 
scales of motion is ambiguous, which is a problem for studies of turbulence: Warn (1986), 
Farge & Sadourny (1989), and Yuan & Hamilton (1994) have all had to invoke the ap- 
proximation of a quadratic energy in their analyses of shallow-water turbulence, which is 
essentially a weak-wave approximation. Second, stability theorems based on the pseudo- 
momentum and pseudoenergy, while easily derived for small-amplitude disturbances (Ripa 
1983), do not appear to extend to finite amplitude (Shepherd 1992). 

What we seek, therefore, is a weak-wave model for shallow-water flow that includes 
the following features: quadratic energy; Hamiltonian structure; appropriate conservation 
laws; quasi-geostrophic or barotropic vortical motion; inertia-gravity waves; and non-trivial 
coupling between the vortical motion and the inertia-gravity waves. 

3. A Hamiltonian weak-wave model 

We begin by transforming the dependent variables in a way that makes the vortical 
and wave motion more transparent in the Hamiltonian structure. In particular, choose 
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q = z ■ (V x v) — frj I H , A = V • v, and rj = h — H as the new dependent variables. Writing 
v = z x VV> + V%, we now non-dimensionalize the variables according to tp = VLip*, x = 



VLx*, q = (V/L)q*, V = (1/L)V*, t = (b/f)U, and rj = HRbq*, where b = B/Vl + B 2 . 
Note that g* = V 2 ^* — brj*, e = Rb, and b < 1. Dropping the asterisks from now on, the 
resulting system can again be cast in the form (2.2) with the state vector a = (q, A,?y) T , 
with J given by 



J = 







I^V(.)) 
^Gfe(-)) 



1+er; ' 

V 2 (-) 



\ 

-V 2 (-) 
> 



(3.1) 



where d(f,g) = f x g y — f y g x , and with the Hamiltonian 

n = J + er/XIVVf + |V X | 2 ) + 2 e?? i • W x V X + 



da;. 



(3.2) 



(Note that in obtaining (3.2), the Casimir invariant f (\gH 2 + gHrj) dx has been removed 
from the original Hamiltonian.) To this point no approximation has been made. 

We now invoke the assumption of a timescale separation, e < 1. To leading order, we 
simply take e = which yields 





(° 








\ 


J = 








-V 2 (- 


) 




Vo 


V 2 (-) 





/ 



n = Jl{m\ 2 + \Vx\ 2 ) + T ^}dx, (3.3a,b) 



with 



-x, 



m 

St] 



V 



1 + B 2 



-bip. 



(3.3c) 



The system (3.3) describes the linearized dynamics 
dq 



dt 



dA - , V 2 T] 



-A, 



(3.4) 
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with no vortical evolution. We now ask: what is the simplest modification that one can 
make to (3.3) to include the nonlinear vortical dynamics, while retaining the Hamiltonian 
structure? Keeping all the 0(e) terms in (3.1) is not an option, because Jacobi's identity 
would then not hold: this is a typical problem in non-canonical Hamiltonian perturbation 
expansions. However, we may now invoke our second scaling assumption, namely that the 
divergence is weak compared with the vorticity: |x| <C \ip\. Note that the first column of 
(3.1) acts on 8H/5q, which to leading order [see (3.3c)] is given by — ip; while the second 
column of (3.1) acts on 5H./5A, which to leading order [see (3.3c)] is given by — x- Thus 
we keep the dominant term in each of the three equations by adding the leading-order part 
of the upper-left entry of (3.1) to (3.3a), namely 



J 



/-ed(q,(-)) \ 

-V 2 (-) 

V o v 2 (-) o / 



(3.5) 



while keeping the quadratic Hamiltonian (3.3b). It may be verified that (3.5) satisfies all 
required properties of cosymplectic forms, including Jacobi's identity. Note, however, that 
although the gravity-wave motion must be weak, it cannot be too weak or the divergence 
equation will be badly approximated by the above procedure. In particular, in the diver- 
gence equation the term arising from the first column of (3.1) must be negligible against 
that arising from the third column; using the fact that the dimensionless gravity-wave 
frequency and q are both O(l), this requires e\ip\ <C The formal scaling regime for the 
validity of this system is therefore: 

e«l, e|V|<|xl<l# (3-6) 

If we absorb the factor e in the nonlinear term by rescaling all the variables by 1/e, then 
(3.3b,c) and (3.5) yields the system 
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This, together with the definition q = V 2 ip — brj, is our weak-wave model. 

The system (3.7) contains both nonlinear vortical dynamics [through (3.7a)] and 
inertia-gravity wave dynamics [through (3.7b,c)]. The coupling between the wave and 
vortical dynamics is seen to come through the linear forcing term bq in the dA/dt equa- 
tion, and through feedback from the wave dynamics onto the dq/dt equation through ijj. 
The system is analogous to the reduced nine-component low-order model of Lorenz (1986, 
Eqn.(2)), which also has all the nonlinearity in the slow equation. 

It is interesting to compare the system (3.7) with the weak-wave model of Farge & 
Sadourny (1989, Eqn.(15)); the latter is given by 

§ = -^9). ^=(lT^")^ + ^-^ V ^)' ^ = -A-d(^). (3.8) 

Evidently (3.8) has additional nonlinear terms, but it is not clear how they arise from 
the original shallow- water system via (3.1). In addition, (3.8) does not appear to be 
Hamiltonian. 

4. Symmetries and conservation laws 

In this section we discuss some of the dynamical implications of the Hamiltonian 
structure of our weak-wave system. The reader may wish to consult Shepherd (1990) for 
further background on some of the points mentioned here. 

The approximation leading to (3.7) consists of keeping the leading order (quadratic) 
approximation to the full Hamiltonian (3.2), together with the simplest approximation to 
the full cosymplectic form (3.1) that incorporates the vortical dynamics yet still satisfies 
the requisite properties, including Jacobi's identity. The resulting weak-wave system (3.7) 
is thereby Hamiltonian by construction. The cosymplectic form (3.5) represents the direct 
product of the non-canonical quasi-geostrophic J, namely —d(q, (•)), and the canonical 
linear- wave J. This kind of structure can be expected to be generic for geophysical fluid 
dynamical systems involving both fast and slow degrees of freedom. 

Hamiltonian structure provides a clear connection between symmetries and conserva- 
tion laws, via Noether's theorem, and allows for the identification of Casimir invariants. 
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The Casimir invariants C are the solutions to J(5C/5a) = 0; in the case of (3.5), this 
condition leads to 

SC , SC SC 

— = function of q, — - = const., — = const., (4.1) 
Sq SA or] 

which implies 

C = J{C{q) + Cl A + c 2 v}dx. (4.2) 

Here C(-) is an arbitrary function and ci, C2 are arbitrary constants. With the exception of 
the second term, which vanishes anyway, the Casimir invariants (4.2) are obvious analogues 
of those present in the original shallow-water system, except that there is no factor of h 
in front of C(q) (this reflects the fact that the advection of q in (3.7a) is solely by the 
non-divergent flow); also, f C(q) dx is identical to the Casimir invariants found in quasi- 
geostrophic dynamics. 

Noether's theorem (see e.g. Benjamin 1984) states that if Sa is an infinitesimal symme- 
try perturbation under which the Hamiltonian is invariant, then a functional M. satisfying 
da = J(5Ai/5a) is a conserved quantity. We consider both zonal translational symmetry 
and rotational symmetry. In the first case, this condition leads to 

-«) = -*• - V V = " A " V JK = - n - (4 ' 3) 

which may be shown to imply (within a Casimir, of course) the zonal momentum invariant 

M = j (yq + v^jdx. (4.4) 

In the case of rotational symmetry with respect to an angle 6>, Noether's theorem leads to 
the angular momentum invariant 

J = J{-\^ + V^)dx. (4.5) 

The terms involving q in (4.4) and (4.5) represent an approximation to the momentum and 
angular momentum of the rotational part of the flow, in vorticity form; they are identical 
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to the momentum and angular momentum invariants of quasi-geostrophic dynamics. The 
terms involving \ in (4.4) and (4.5) represent the momentum and angular momentum of 
the divergent part of the flow. 

5. Pseudoenergy, pseudomomentum, and nonlinear stability 

A general property of Hamiltonian systems is that any basic state possessing a sym- 
metry invariance is a constrained extremal of the integral invariant corresponding to that 
symmetry, the constraint being imposed through a suitable Casimir. This means that one 
may construct exact, finite-amplitude invariants that are quadratic to leading order in the 
disturbance to such basic states. For example, with steady basic states one can combine 
the energy and Casimir invariants to obtain a so-called pseudoenergy invariant. In the 
same way, with zonally symmetric basic states one can combine the zonal momentum and 
Casimir invariants to obtain a pseudomomentum invariant. Such invariants can be used 
to derive nonlinear stability theorems. See Shepherd (1990, 1994) for further background. 

(a) Zonally symmetric basic states 

We first consider disturbances {<f , A', rj'} to a steady, zonally symmetric basic state 
{qo(y), A (y), r)o(y)}. In this case we can use the energy (3.3b), Casimir (4.2), and zonal 
momentum (4.4) invariants to construct the combined invariant 

A = (H + C + aM) [q + q', A + A', r) + ri']- (H + C + aM) [go, A , 770], (5.1) 

where a is kept free. The Casimir C is determined by the extremal condition 8A = 0, 
which leads to 

C'(qo) = V>o - ay, A = 0, 770 = 6(1 + B 2 )^ . (5.2a,b,c) 

(Note that c\ and C2 can be set to zero without loss of generality by imposing a suitable 
gauge condition on the potentials ipo an d Xo-) Conditions (5.2b,c) express geostrophic 
balance, which is required for a steady solution. Condition (5.2a) determines the functional 
dependence of C(-); thus the right-hand side of (5.2a) must be interpreted as a function of 
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qo. If we introduce the function G a (-) denned by G a (qo(y)) = ipo(y) — ay, then A can be 
written 



rqo+q q / ^1 

-4>o<f + J G a (Z)dZ + ayq' + arj' j dx 

^ + (f ) 2 + ( v^w + -v^^) 2 + [i - o-d + *>] (f) 2 ) 

+ ^ [<?a(?o + 0- (5-3) 

This expression is seen to be positive definite if 

(U o + a)^-<0 Vy and a 2 < 1+ 1 g2 , (5.4a,b) 

using dGa/dqo = [dG a / dy]/ (d%/ dy) = — {Uq + a) / (dqo / dy) , where £7o = —dipo/dy. (Note 
that under the condition (5.4a), G a (-) is guaranteed to be single- valued.) Since A is an 
exact invariant of the nonlinear dynamics, it follows using standard methods that whenever 
a can be chosen such that (5.4a,b) are satisfied, the basic state is nonlinearly stable. In 
particular, a normed stability theorem can be established!, which holds for disturbances 
of arbitrarily large magnitude. 

It was Arnol'd (1966) who first introduced this approach [building on earlier work of 
Fj0rtoft (1950)] and applied it to the barotropic system. For that system, (5.4a) alone is a 
sufficient condition for nonlinear stability. The same result applies to the quasi-geostrophic 
system. Our stability theorem is thus seen to be the analogue of Arnol'd's theorem for the 
weak-wave model, with the additional condition (5.4b) arising from the model's gravity 
wave (unbalanced) dynamics. 

Ripa (1983) followed the above methodology for the full shallow-water system, and 
showed that the second variation of the corresponding invariant A was positive definite if 
an a could be found such that (in terms of our present variables) 

(U o + a)^<0 Vy and ( a + U ) 2 <±±^ Vy, (5.5a,b) 
dy 1 + B z 



f Technically this requires that the left-hand side of (5.4a) be bounded away from zero 
and minus infinity. 
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where q = (b+V 2 ipo) / (1+770) is a scaled potential vorticity. Note that under the rescaling 
implicit in (3.7), 1 770 1 = 0(e) <C 1 so q b + q and the right-hand side of (5.5b) reduces 
to the right-hand side of (5.4b). Thus the stability theorem resulting from conservation 
of A in the weak-wave model is analogous to that in the full shallow-water model. The 
key difference is that the higher-order terms in the shallow-water form of A cannot be 
controlled (it appears that only the second variation of A can be shown to be positive 
definite), so condition (5.5) is only a statement of small-amplitude (linear) stability. 

The fact that the weak-wave stability theorem is fully nonlinear means that it can 
be used to obtain rigorous upper bounds on the nonlinear saturation of instabilities in 
this model, following the method of Shepherd (1988). In Section 6 it is shown how this 
technique places rigorous bounds (within the context of this weak-wave model) on the 
emission of gravity waves from unstable parallel vortical flows. 

(b) Axisymmetric basic states 

We now consider disturbances to a steady axisymmetric basic state {qo (r), A (r), i] (r)} 
where x = rcos6>, y = rsin6>. In this case we proceed exactly as above, but using the an- 
gular momentum invariant (4.5) in place of the zonal momentum. Defining V according 
to (5.1) but with M. replaced by J, the extremal condition 5V = leads to 

2 

C'(q ) =^o + «y, Ao = 0, rio = 6(1 + B 2 )^ - (5.6a,b,c) 

Introducing the function G a (-) defined by G a (<?o( r )) = V ; o( r ) + ar 2 /2, it follows that V 
can be written 



T= M^|VV>? + |Vxf + :^)-r l 



+ 



rqo+g r 2 q / ^ 

W + J G o (0 dti - a-q' + ctr/-±j dx 

G a (qo + Z)-G a (q j\dt}dx. (5.7) 
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As with A, V is an exact invariant of the nonlinear dynamics. The expression (5.7) is seen 
to be positive definite if 

(U + ar)^>0 Vr and a 2 < 1 -- , (5.8a,b) 
dr r z (l + B z ) 

where Uq = difto/dr. However, it turns out that for physically realizable flows in an 
unbounded domain, (5.8a,b) can never be satisfied! To see this, first note that taking 
r — > oo in (5.8b) implies that a = 0. Using (5.6c), q = d 2 ip /dr 2 + (l/r)dip /dr — B 2 ip . 
Therefore (5.8a) requires 

7 dr dr 7o ^ <^ r <^ r3 dr dr\r dr J V dr / J 



using the conditions dllo/dr — > as r — > oo and £/o — > as r — > 0. But 



<ir <ir Vr dr / ./n ^ r dr dr 2 r 2 V dr 



J_ /d^ \ 2 i 
r 2 V dr / J 

= _r#od/l*. r _ Tlf^fdr, (5.10) 
7 dr dr Vr dr / 7 r 2 V dr / 

since £/o = O(r) as r — > for finite vorticity. Hence the second integral on the right-hand 
side of (5.9) is negative definite, which implies that (5.8) can never be satisfied. 

In the barotropic or quasi-geostrophic case, (5.8a) alone is a sufficient condition for 
nonlinear stability, with no restriction on a. In the shallow-water case, even a small- 
amplitude stability theorem cannot be obtained by this approach for unbounded domains 
(see Appendix). So we conclude that the weak- wave model captures the essential physics 
of the full shallow-water dynamics with regard to stability theorems of the Arnol'd type. 

6. Rigorous upper bounds on gravity- wave emission from unstable parallel 
shear flows 



Suppose we are given a parallel shear flow {qi(y), 7/i(y)} which is known to be unsta- 
ble. In general we expect the instability, even if it is vortical in nature, to excite gravity 
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waves. For example, Ford (19946) has studied shear-layer ( Ray leigh- type) instability in the 
shallow-water equations, and has found the generation of gravity waves to be ubiquitous. 
An important question is whether there are any a priori upper bounds on the amplitude 
of such gravity-wave emission. 

Such bounds can be obtained using our nonlinear stability theorem. Following Shep- 
herd (1988), we may regard the initial unstable flow as a finite-amplitude disturbance to 
some nearby stable flow, and then invoke the nonlinear stability theorem. Since steady 
flows are necessarily balanced (with no gravity- wave component), any bound on the dis- 
turbance will also provide a bound on the gravity-wave activity. We provide an explicit 
demonstration of this by deriving a bound on the kinetic energy of the divergent flow, 



f\\V X \ 2 dx. (6.1) 



Let us assume, for the sake of simplicity, that the unstable flow is perturbed infinitesimally 
at t = 0. (It would be straightforward to include the effect of a finite initial perturbation.) 
Now, introduce a basic state {qo(y), Voiv)} that is stable by criteria (5.4). At t = 0, the 
disturbance {<f , A', r]'} relative to this basic state is given by r]'(0) r\\ — i]o : etc., and the 
invariant A is approximately given by 

For t > 0, we note that x = x' (since the basic state has no divergent component); then 
from (6.1) and (5.3), together with the conservation of A in time, we obtain the inequality 

f m < - -4(0) 

td[ V " 1 - «2(1 + S2) " 1 _ «2 (1 + B 2 ) > ^ 

into which (6.2) may be substituted. (Note that the factor 1 — a 2 (l + B 2 ) is always positive, 
by virtue of (5.4b).) In this way we obtain an a priori upper bound on the divergent kinetic 
energy Sd(t) for all time, which is a functional of the initial unstable flow and of the chosen 
stable basic state (including the choice of a) . This places a rigorous bound on the emission 
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of gravity waves from the unstable flow. Since the basic state is arbitrary, one may try 
to minimize the bound over all possible choices of stable basic state (cf. Shepherd 1988). 
In the context of the present weak-wave model, this analysis has been performed by Nore 
(1994) for the Bickley jet and the tanh(y) shear layer. 

For an initial flow that is marginally unstable, the above procedure yields a bound 
that is a function of the supercriticality — implying weak gravity-wave emission. Argu- 
ments from statistical mechanics (Warn 1986) would suggest, to the contrary, that once a 
steady flow goes unstable, its energy will be partitioned roughly equally into both vortical 
and gravity-wave components. The nonlinear stability theorem therefore provides definite 
limits on the extent of such "thermalization" in the weak- wave model (cf. Shepherd 1987). 

7. Linear eigenvalues and instability 

The shallow- water equations allow mixed vortical/gravity- wave instabilities that are 
filtered by the quasi-geostrophic equations. An important test of our weak-wave model is 
the extent to which it can capture such instabilities. As a preliminary step, we consider 
the instability properties of two simple basic flows. 

(a) Zonally symmetric basic state 

We first consider the case of a piecewise-constant zonally symmetric potential vorticity 
distribution with a single jump at y = 0. As such a flow has no intrinsic length scale, we 
may take B = 1 without loss of generality, and introduce the basic state 

(2Q, y>0 (-2Q + Qe~y, y>0 

Qo = \ V>o = < , (7-1) 

l0, y<0 ( -Qe*, y<0 

with 7]o = y/2ipo. Again without loss of generality, we take Q > 0. The velocity field 
associated with this basic state is given by Uo(y) = Qe _ ' y ', namely an eastward symmetric 
jet with a maximum at y = 0. 

Physically this problem consists of gravity waves in the upper and lower half-planes, 
coupled to a potential-vorticity disturbance at y = 0. In each of the two half-planes we 
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introduce a normal-mode form for if;', proportional to e l { kx + e -y-^ t ) ^ anc [ insist on evanes- 
cence as \y\ — > oo. The frequency u must of course be the same in the two half-planes; 
the same condition applies to the (real) zonal wavenumber k, by virtue of the kinematic 
condition that v' = ip x must be continuous across y = 0. Away from y = 0, we have q' = 
and so (3.7b,c) — with B = 1 — yield 

u 2 = ^(l + k 2 + £ 2 ), (7.2) 

which is just the dispersion relation for inertia-gravity waves. We obtain a dynamical 
matching condition at y = by integrating the linearized potential vorticity equation 
(3.7a), namely q' t + Uoq' x + ip' x qo y = 0, across y = 0. This yields (ui — Uo(0)k)[£] = —ik[qo], 
where [•] denotes the jump across y = 0. Evidently a jump in t is required to balance the 
jump in qo. Since (7.2) implies that I 2 must be the same on either side of y = 0, this means 
that i must change sign. Letting I denote the value for y > 0, the dynamical matching 
condition may be written 

(oj - Qk)£ = -iQk. (7.3) 

Equations (7.2) and (7.3) define the linear disturbance problem. Note that ii > for 
evanescent disturbances; it may be shown that this requires the intrinsic frequency Cj = 
oo — Qk to have a negative real part. For quasi-geostrophic dynamics, the condition (7.2) 
is replaced by the condition £ 2 = —k 2 — 1, which when substituted into (7.3) leads to the 
single (stable) solution branch 

^= Qk ~^w- (7 ' 4) 

Condition (5.4a) is satisfied (in the appropriate limiting sense at y = 0) if a < —Q, 
so the basic state (7.1) is provably stable by our criterion (5.4) for Q < Q c = l/y/2. This 
stability condition has a simple physical interpretation. From (7.2), the phase speed of the 
gravity-wave mode is always greater than l/v^, but approaches that value in the limit 
k — > oo. From (7.4), the phase speed of the quasi-geostrophic mode is always less than Q, 
but approaches that value in the limit k — > oo. So when Q < Q c , the phase speeds of the 
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two modes are always separate; but when Q > Q c , they cross for sufficiently large k. Our 
stability condition (5.4) is therefore a condition for the absence of a mode crossing. 

In terms of dimensional variables, the critical value corresponds to [q] = 2f. This lies 
outside the formal domain of validity of the weak-wave model, since it requires order-unity 
relative changes in the layer depth. 

Combining (7.2) and (7.3) leads to the quartic equation 

2u 4 + 4Qku 3 + ([2Q 2 - l]k 2 - l)cD 2 + Q 2 k 2 = 0. (7.5) 

We can analyze (7.5) asymptotically in the limits of small and large wavenumber k. Note 
that the coefficient of Cj 2 changes character depending on the value of Q: for Q < Q c , 
it is always negative; for Q > Q c , it changes sign at some value of k. This latter case 
corresponds, as anticipated, to a mode crossing, as shown below. 

For k <C 1, the method of dominant balance leads to the two stable solutions 

2u> 4 ~ Cj 2 =>- Cj ~ — , u ~ ■= (7.6) 

V2 ^2 



and 



uj 2 ~Q 2 k 2 u~-Qk, (7.7) 



In each case there are two real roots but we have kept only the one with Qj < 0. The root 
(7.6) evidently corresponds to the small-/c limit of a pure gravity-wave mode (7.2), while 
the root (7.7) corresponds to the small-/c limit of the quasi-geostrophic mode (7.4). (The 
expression for u in (7.7) obtains on expanding (7.5) in k.) 

For k ^> 1, the method of dominant balance leads to the solutions 

2u 4 + AQku; 3 ~ (1 -2Q 2 )k 2 u 2 =>- cu~-QA;±-^, u ~ ±-^= (7.8) 

v 2 v 2 



and 



(1 - 2Q 2 )k 2 u 2 ~ Q 2 k 2 => u ~ ±-—2==, ur^Qk. (7.9) 

V 1 - 2Q 
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There are two distinct cases here. When Q < Q c , then only the negative root of (7.8) 
has Oj < 0; it is evidently the large-/c limit of a pure gravity-wave mode (7.2). In this 
case another stable solution is provided by the negative root of (7.9); it reduces to the 
large-/c limit of the quasi-geostrophic mode (7.4) for Q <C 1, and generally has the same 
characteristic of the intrinsic frequency asymptoting to a constant. When Q > Q c , on the 
other hand, (7.9) implies an unstable solution with uji ^ 0; in this case both roots of (7.8) 
have Cj < 0. 

Therefore, for Q < Q c , when our stability theorem predicts stability for this basic flow, 
we find that both quasi-geostrophic and gravity-wave neutral modes can be identified, and 
are connected smoothly (i.e. arise from the same dominant balance) in the small-/c and 
large-/c limits: there is no mode crossing. But for Q > Q C: there is a dramatic change in 
the character of one of the solutions: the small-/c quasi-geostrophic mode metamorphoses 
into a gravity-wave mode in the large-/c limit, and an unstable solution appears. 

Ford (1993) has studied the analogous basic state (with piecewise-constant q rather 
than piecewise-constant q) in the context of the shallow- water equations, obtaining distur- 
bance equations analogous to (7.2)-(7.3) but with the full basic-state height field entering 
the problem and with a Doppler shifting of the gravity-wave frequency. (Both of these 
effects are neglected in the weak-wave model.) Ford found a similar transition in the 
character of the neutral modes at a critical value of Q, corresponding to the critical value 
predicted by Ripa's stability criterion (5.5); for Q larger than the critical value, he found 
mode crossing and instability. In terms of our non-dimensionalization, Ford's critical value 
of Q (based on the jump in relative vorticity) is 3/ (4\/2), which is somewhat less than our 
Q c . Moreover, in the shallow-water case the growth rates are exponentially small for large 
k, whereas here they are of order unity. Both these discrepancies reflect the fact that the 
instability occurs outside the range of validity of the weak-wave model. 
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(b) Axisymmetric basic state 



We now consider the case of a piecewise-constant axisymmetric potential vorticity 
distribution with a single jump at r = 1: 

(Q, r<l ( -Q + BQK 1 (B)I (Br), r<l 

% =\ , B 2 %=\ , (7.10) 

U, r>l ( -BQIi(B)K (Br), r>l 

with rjo = (1 + B 2 )bipQ. Here I m and K m are modified Bessel functions of the first and 
second kind, of order m. The velocity field associated with this vortex is continuous in r, 
and has a maximum at r = 1 given by Uq{1) = QK\{B)I\{B). 

As in the zonally symmetric case, this problem consists of gravity waves inside and 
outside the vortex, coupled to a potential vorticity disturbance at r = 1. Away from 
r = 1, we introduce a normal-mode form for ip', proportional to f^^i™-®-^ 1 ) where 
f(r) = Im{ir) for r < 1 and f(r) = I m {l)K m (^r) / 'K m (^) for r > 1 (so f(r) is continuous 
across r = 1). Taking q' = for r ^ 1, (3.7b,c) yield the inertia-gravity wave dispersion 
relation 

(1 + B 2 )uj 2 = B 2 - 7 2 . (7.11) 

As above, we obtain a dynamical matching condition at r = 1 by integrating the linearized 
potential vorticity equation across r = 1. After some algebra, this can be shown to reduce 
to 

to = QmlK^hiB) - K m ( 7 )/ m ( 7 )]. (7.12) 

Equations (7.11) and (7.12) define the linear disturbance problem. The quasi-geostrophic 
problem consists of (7.12) alone with 7 replaced by B, which gives a single (stable) mode. 
In the limit B — > 0, this solution reduces to Kelvin's (1880) result to = Q(m — l)/2. 

The absence of a stability theorem in the axisymmetric case suggests the existence 
of mode crossing between the gravity-wave and quasi-geostrophic modes for all parameter 
values. That this is indeed the case can be seen as follows. In the limit m ^> max{S, 1}, 
(7.12) suggests that uj = O(m); then (7.11) implies 7 ~ ±iy/l + B 2 uj. Since K m (^)I m (^) = 
0(1/7) for 7 » 1, the second term of (7.12) is much smaller than the first, which is 
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consistent with the assumption that u = 0(m). Substituting this expression for 7 back 
into (7.12), an imaginary component to u is implied: in particular, 



u r ^Q m K 1 (B)I 1 (B), u t ~ ±|Qto \j m Q 2 m 2 K 2 {B)I 2 {B){l + B 2 ) - B 2 

(7.13) 

where J m is the Bessel function of the first kind, of order to. We thus find instability for 
sufficiently large to, for all parameter values. In the limit of small Q and large to, we may 
make the substitution sech/3 = QKi(B)Ii(B)\ / l + B 2 in (7.13), and use the asymptotic 
approximation 

e m(tanh/3-/3) 

J m (m sech/3) ~ - (7.14) 
V27tto tanhp 

Together with the identity e 2/3 = (1 + tanh/?)/(l — tanh/3), this yields an approximate 
growth rate 



Ui ~ — 



q e 2mtanh/3 /l_ t anh/3\ m 



4 tanh /3 V 1 + tanh (3 
Now using the small-Q expansion 



J • (7-15) 



tanh/3 = yj 1 - sech 2 /? ~ 1 - - Q 2 Kf(B)I 2 (B)(l + B 2 ) + . . . , (7.16) 



we obtain 



I e 2m (1 QKi(B)/i(5)(l + S 2 ) 1 / 2 ) 2 ". (7.17) 



The derivation of this expression, valid for small Q and large to, was provided by R. 
Ford (personal communication, 1996). The Froude number based on the velocity at the 
vortex edge is F = (1 + S 2 ) 1 / 2 t/" (l), namely F = QK 1 {B)I l {B){l + B 2 ) 1 / 2 , so we see 
from (7.17) that the growth rate scales as F 2m . This is in accord with the instability 
analysis of Ford (1994a) for the analogous basic state in the shallow-water equations (with 
piecewise-constant q rather than piecewise-constant q); Ford obtained 

"•~i e2<m - i, (?r < 7 - i8 > 
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in the limit of large m and small F. The factor of Q difference between (7.17) and (7.18) is 
accounted for by the different nondimensionalization of timescales, and the factor of 2 2m 
difference by the fact that Ford's vortex-edge velocity is 1/2 rather than 1 (so that there 
is a factor of two difference in the definition of F). The only remaining difference between 
(7.17) and (7.18) is a factor of e 2 , the origin of which is unclear. 

When B < 1, then + B 2 ) 1 ' 2 -> 1/2 and the approximation (7.15) is 

still valid for Q < 2 (in the large-m limit). In this case, the weak- wave model predicts 
a growth rate of 0(1) when Q = 0(1). This is very much larger than the shallow-water 
growth rate of 0(F 2m ) (see (7.18)) for B <C 1 and Q = 0(1), and indicates that in this 
regime the weak- wave model is outside its range of validity. The problem is that although 
B <C 1 for the vortex, the instability involves a gravity wave with a very long wavelength, 
whose frequency is not much larger than the rotation rate of the vortex. This illustrates 
the pitfalls of a formal scaling analysis based on a single length scale. 

8. Conclusion 

Atmospheric and oceanic flows can often be characterized as consisting of slow nonlin- 
ear vortical motion evolving in the presence of fast but weak linear-wave motion. We seek 
a reduced dynamical model that can describe this situation, and that has the following 
features: quadratic energy; Hamiltonian structure; appropriate conservation laws; quasi- 
geostrophic or barotropic vortical motion; fast waves; and non-trivial coupling between 
the vortical motion and the fast waves. In this paper we have derived such a "weak-wave" 
model in the context of the /-plane shallow-water equations, where the fast waves are 
inertia-gravity waves. 

Although our model is appealing from the Hamiltonian perspective, it has at least 
two shortcomings. The first is that our method of approximation is not unique: if we 
had chosen different prognostic (dynamical) variables, we would have obtained a some- 
what different system. This reflects the fact that our model is not based on a rigorous 
perturbation expansion. Against this we can only say that the problem of combining rig- 
orous perturbation expansions with non-canonical Hamiltonian structure remains an open 
problem, which nobody has yet been able to solve. 
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The second shortcoming of our model is that, for formal self-consistency, the fast waves 
must not be too weak. This might cast doubt on the utility of the model for studying 
mixed vortical/gravity- wave instabilities, and the breakdown of balanced dynamics. To 
determine whether the model is deficient in these respects, we have analyzed its stability 
and instability properties. 

With regard to stability, the model possesses an Arnol'd-type stability theorem for 
zonally symmetric vortical basic flows, (5.4), that corresponds to the analogous theorem 
(5.5) for shallow-water dynamics, and that has the same physical interpretation: a basic 
flow is stable if its velocity is everywhere subcritical with respect to the gravity- wave phase 
speed. However, in our case the stability theorem is fully nonlinear. This fact permits the 
establishment of rigorous upper bounds on gravity-wave emission from unstable parallel 
shear flows. As is the case with shallow- water dynamics (see Appendix), an Arnol'd-type 
stability theorem cannot be obtained for axisymmetric vortical basic flows. This is in 
striking contrast with the quasi-geostrophic case — where any vortex with a monotonic 
potential-vorticity profile is nonlinearly stable - - and reflects the destabilizing role of 
gravity waves in our model. 

With regard to instability, we have examined the linear dynamics of disturbances to 
two piecewise-constant potential-vorticity distributions. In the case of the zonally symmet- 
ric basic state, corresponding to an exponential jet, the nature of the eigenvalue problem 
depends on the jet strength. When the jet is subcritical, the modes are everywhere neutral 
and can be clearly identified with quasi-geostrophic and gravity-wave modes. When the 
jet is supercritical, however, the eigenvalues collide (mode crossing) and instability occurs 
for sufficiently short zonal wavelengths. This behaviour corresponds qualitatively to that 
seen in the shallow-water system (Ford, 1993). In the case of the axisymmetric vortex, 
instability is found for all parameter values: gravity waves destabilize a vortex that would 
be stable under balanced dynamics. Once again, this is qualitatively consistent with the 
behaviour of the shallow- water system (Ford, 1994a). It is also consistent with the study 
of Broadbent & Moore (1979), who found a destabilization of axisymmetric vortices by 
coupling to acoustic waves in the context of the two-dimensional compressible equations. 
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The gravity-wave equation implied by (3.7) is 



d 2 A 




B 2 - V 2 
1 + B 2 




(8.1) 



dt 2 



Although the zeroth moment of the right-hand side of (8.1) clearly vanishes (for a compact 
vorticity distribution), the first moment does not appear to vanish. This is in contrast 
to the shallow-water system, which has vanishing zeroth and first moments (Ford et al. 
1996), and suggests that our weak- wave model might overestimate the forcing of gravity 
waves by vortical motions (Ford, personal communication, 1996). However, the d(ift, V 2 ip) 
term in (8.1) can be put into an explicitly quadrupolar form, which leaves only the d(^, rj) 
term as a dipolar source. For weak waves, r\ is approximately given by 6(1 + B 2 )ip (this is 
the quasi-geostrophic approximation), so this term may be expected to vanish to leading 
order in e. Thus it may well be that the gravity-wave source term in (8.1) has reasonable 
scaling properties in e. Further investigation of this point, through numerical simulations, 
would be of considerable interest. 

Because the energy of our weak-wave model is quadratic, and the potential vorticity 
is a linear function of the other variables, this means that one can apply standard spectral 
arguments to deduce statistical mechanical equilibria and turbulent spectral regimes. In 
fact, Warn (1986) and Farge & Sadourny (1989) have already done so. Although they were 
studying the shallow-water equations, they made the approximation of a quadratic energy 
and a linear potential vorticity in order to facilitate their analysis. It follows that their 
results apply directly to our weak-wave model. 
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Appendix. On the non-existence of Arnol'd-stable circular 
vortices in unbounded shallow-water flow 

Using the Arnol'd-stability methodology followed in Section 5, Ripa (1987) derived a 
stability condition for axisymmetric basic flows in the context of the shallow-water equa- 
tions. In terms of dimensional variables, Ripa (1987) showed that the small-amplitude 
(quadratic) approximation to the combined energy/angular-momentum/Casimir invariant 
V was positive definite if there existed a constant a such that 

(V + ar) > Vr and (V + arf < gH Vr. (A.la,b) 

Here V(r) is the basic-state azimuthal velocity, Q(r) = (/ + dV/dr + V/r)/H is the basic- 
state potential vorticity, and H(r) is the basic-state total height, satisfying the gradient- 
wind balance condition gdH/dr = fV + V 2 /r. Conditions (A.la,b) are clearly analogous 
to the stability conditions (5.8a,b) for our weak-wave model, after the appropriate non- 
dimensionalization. 

However, it turns out that for physically realizable flows in an unbounded domain 
(A.la,b), like (5.8a,b), can never be satisfied. (Ford (1994a) has provided a heuristic 
argument to this effect for the special case of monotonic potential-vorticity profiles.) We 
confine our attention, as in Section 56, to localized vortices with dV/dr — > as r — > oo. In 
that case H(r) = o(r 2 ) as r — > oo, and so the only way that (A. lb) can be satisfied in an 
unbounded domain is by taking a = 0. In that case the angular momentum invariant is 
not being used, and as Ripa (1987) argues, Arnol'd stability therefore cannot be proven. 
The argument involves an appeal to the result of Andrews (1984): if V consists only of the 
energy and a Casimir invariant, neither of which are altered by a spatial translation of the 
vortex, then the basic state cannot be a true extremum of V for arbitrary perturbations, 
and V cannot be sign definite. 
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There are, however, some subtleties involved with the use of Andrews' theorem in an 
unbounded domain (Carnevale & Shepherd 1990), and it is therefore desirable to establish 
the non-existence of Arnol'd-stable circular vortices by direct methods, as was done in 
Section 56 for the weak-wave model. By definition, rV(r) = f Q (QH — f)r' dr'; this implies 
that V = 0(r) as r — > (for finite vorticity) and that dV/dr — > as r — > oo (for a localized 
vortex, with QH — f — > as r — > oo). Taking a = as argued above, we now try to show 
that (A.la,b) can never be satisfied. Noting that 




(A.3) 

Now, (A. lb) requires V 2 < gH, in which case the right-hand side of (A.3) is negative 
definite; but (A. la) requires the left-hand side of (A.3) to be positive definite. It follows 
that Ripa's stability conditions (A.la,b) can never be satisfied in an unbounded domain. 
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